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The study of high-index dielectric nanoparticles currently attracts a lot of attention. They do 
not suffer from absorption but promise to provide control on the properties of light comparable to 
plasmonic nanoparticles. To further advance the field, it is important to identify versatile dielectric 
nanoparticles with unconventional properties. Here, we show that breaking the symmetry of an all¬ 
dielectric nanoparticle leads to a geometrically tunable magneto-electric coupling, i.e. an omega-type 
bianisotropy. The suggested nanoparticle exhibits different backscatterings and, as an interesting 
consequence, different optical scattering forces for opposite illumination directions. An array of 
such nanoparticles provides different reflection phases when illuminated from opposite directions. 

With a proper geometrical tuning, this bianisotropic nanoparticle is capable of providing a 27r phase 
change in the reflection spectrum while possessing a rather large and constant amplitude. This 
allows creating reflectarrays with near-perfect transmission out of the resonance band due to the 
absence of an usually employed metallic screen. 

PACS numbers: 42.25.-p, 78.67.Bf, 78.20.Bh, 45.20.da,42.25.Fx 

Nanophotonics has attracted enormous research inter¬ 
ests due to its potential to control light-matter interac¬ 
tion at the nanoscal^il®. Being usually connected with 
a strong light confinement in metal-dielectric and plas¬ 
monic structures, nanophotonics offers remarkable op¬ 
portunities due to the local field enhancementP^^. How¬ 
ever, applications of plasmonic nanophotonics are often 
limited by Ohmic losses at optical frequencieJ^H^. An al¬ 
ternative strategy in nanophotonics is to use high-index 
dielectric materials for building blocks that control light- 
mater interactiorP^®^. Dielectric nanoparticles have at¬ 
tracted a considerable attention due to their appealing 
application^I^H^. 

For these nanoparticles to be versatile, usually the mul¬ 
tipolar composition of their scattered field is tailored. 

A high-index dielectric nanoparticle that supports both 
electric and magnetic resonant responses where the ex¬ 
act composition of the two contributions can be tuned 
by means of geometrical modifications has been recently 
suggestecP^H^. Interesting optical features such as direc¬ 
tional scattering pattern with zero backscattering or an 
optical pulling force can be achieved by proper tuning 
of these response^^^^^. Nanoparticles with higher multi¬ 
polar responses, e.g. electric and magnetic quadrupoles 
may provi de an even more control on the scattering prop- 
ertie^SM] compared to nanoparticles with only 

dipolar responses. However, they require more compli¬ 
cated theoretical considerations and a range of dissim¬ 
ilar material^^. Therefore, it is a challenge to bring 
more control on the scattering features of nanoparti¬ 
cles without adding such complexities. Here, we use 
hianisotropy to control the scattering properties of high- 
index dielectric nanoparticles in the context of dipolar 



FIG. 1: (a) Schematic of the proposed bianisotropic high- 
index dielectric nanoparticle, (b) A cut in the yz-plane of the 
proposed nanonatenna. The dielectric cylinder has a diameter 
D = 300 nm and height H — nm with a cylindrical 
hole inside it that has a diameter Do = 150 nm and height 
Hq — 10 — 140 nm. We assumed that the refractive index of 
the material from which the cylinder is made is n — 3.5. The 
ambient material shall be air with a refractive index of unity. 


response d^^ l ^^ i Note that we are only interested in the re¬ 
ciprocal omega-type bianisotropy and all the other types 
of bia nisotr opy (i.e. chiral, Tellegen, and moving) are 
absenfP^J^. We show that an omega-type bianisotropic 
response can be achieved by breaking the symmetry of 
a cylindrical nanoparticle, i.e. a high-index dielectric 
cylinder with partially drilled cylindrical air hole [see 
Figs.j^a) and [^b))]. 

A reciprocal omega-type bianisotropic nanoparticle 
supports electric and magnetic dipole moments with dif¬ 
ferent strength when illuminated by a plane wave in for¬ 
ward or backward directions [see Fig.[^a)]. In particular, 
the bianisotropy causes significant difference in backseat- 
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FIG. 2: (a) and (c) Total scattering cross sections and 
contributions from different Cartesian multipole moments as 
a function of frequency for forward (+) and backward (—) il¬ 
lumination directions; electric dipole moment (red dashed 
line), magnetic dipole moment (blue dashed line), elec¬ 
tric quadrupole moment (green dashed line) and magnetic 
quadrupole moment QM^ (purple dashed line). The height 
of the air hole is Hq = 100 nm. (b) and (d) normalized mag¬ 
netic field distributions at 365 THz [gray dashed line in (a) 
and (c)] for both opposite illumination directions. 


tering responses, and consequently in the optical scat¬ 
tering forces exerted on the proposed nanoparticle for 
opposite illumination directions. The observed magneto¬ 
electric coupling can be tuned by controlling geometric 
parameters of the nanoparticle. We demonstrate that an 
infinite periodic planar array of such lossless nanoparti¬ 
cles (a metasurface), provides different resonant reflec¬ 
tion phases for different illumination directions. This is 
due to the proposed magneto-electric coupling, and it is 
certainly not possible with only an electric and/or mag¬ 
netic response. We also demonstrate that the investi¬ 
gated nanoparticle is an excellent candidate for unit-cells 
of reflect arrays. Indeed, this nanoparticle in its reflec¬ 
tion resonance band, provides a 27r phase change for a 
proper geometrically tuning while it maintains a con¬ 
siderable large and constant amplitude, similar to Huy¬ 
gens’ metasurface^^. Moreover, since there is no metal¬ 
lic ground plate in the proposed reflectarray, it will be 
transparent out of its resonance band. 


I. BIANISOTROPIC NANOPARTICLES 


The geometry of the proposed bianisotropic high-index 
dielectric nanoparticle is shown in Figs, [^a) and Bb). 
The suggested nanoparticle is not symmetric with re¬ 
spect to the forward and backward illumination direc¬ 


tions [see Fig. [^a)]. Material and geometric parame¬ 
ters of the nanoparticle are given in Fi g. fTj Using the 
multipole expansion of the scattered flelcP^, we compute 
the total scattering cross section and the contribution 
of each multipole moment when the nanoparticle is illu¬ 
minated by a plane wave in the forward and backward 
directions [see Figs.j^a) andj^c)]. The higher order mul¬ 
tipole moments, i.e. electric quadrupole Q and magnetic 
quadrupole QM moments are negligible below 450 THz 
[the frequency range above 450 THz, signifying the upper 
frequency for the applicability of the dipole approxima¬ 
tion, is marked by the gray shadowed areas in Figs. da) 
and[^c)]. In the following, we only concentrate on this 
frequency range (i.e. lower than 450 THz where only elec¬ 
tric and magnetic dipole moments are dominant). The 
corresponding electric and magnetic dipole mo¬ 
ments of the nanoparticle for both illumination directions 
[see Fig. are given by: 

^ ( 1 ) 

eo 

Zom^ = ± ( 2 ) 

where eo is the free space permittivity and Zq is the char¬ 
acteristic impedance of free space, o^ee, <^em, <^mm 
are individual polarizabilities of the nanoparticle. They 
represent electric, magneto-electric, electro-magnetic, 
and magnetic couplings in the proposed nanoparticle, re¬ 
spectively. That is, the interactions b etwee n the electric 
and magnetic fields and polarization d^^FS I Notice that 
in order to distinguish the illumination direction of the 
nanoparticle, we use the ± sign where the plus/minus 
sign corresponds to the propagation direction of the in¬ 
cident plane waves in forward/backward direction (the 
time dependency is assumed to be In order to ex¬ 

plain the underlying physical mechanisms of the scatter¬ 
ing response of the bianisotropic nanoparticle, we start 
with the definition of extinction cross section for 

a bianisotropic nanoparticle in dipole approximation for 
both illumination directions, i.e. FF^FO ] 


_ ext 

^ext T 

J-0 

-|Re (Einc X X ■ nds 

lo 

— klm. (o^ee i i T <^mm) 5 (3) 

where Iq = |Einc|/2^o is the intensity of the incident 
plane wave. |Einc| and are the amplitude of the 

incident and scattered electric fields for both illumina¬ 
tion directions (±). is the extracted power by 

the nanoparticle, k = cj/c is the wavenumber for an 
angular frequency uj. For a reciprocal nanoparticl^^, 
i.e. a^rn = —<ame: wc cau conclude that the extinc¬ 
tion cross section for both illuminations are identical, i.e. 
Cext = Cit = Cext and is given by 

^ext kJ-IH (ttee n^mm) 5 


( 4 ) 
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On the other hand, we know for the proposed nanopar¬ 
ticle that the scatterin g cross sections for both illumina¬ 
tions can be written 


= 


-Psca 

lo 


I Re X ■ nds 

lo 


— ^l^ee ^ ^em| "t“ |^me ^ ^mm 


), ( 5 ) 


where, is the radiated or scattered power by the 
nanoparticle. The extinction cross section Cgxt for a 
lossless reciprocal nanoparticle is identical to the scat¬ 
tering cross section Cgca? he. Cext = C^ca due to the 
fact that the absorption cross section is zero Cabs = 
C'ext ~ Csca = 0- This explains why the scattering 
cross sections are identical for the proposed nanoparticle 
when illuminated from opposite directions, i.e. Cgca = 
C's^a — ^Pca black solid lines in Figs. a) and[^c)]. 
Notice, the equality between extinction and scattering 
cross sections for the lossless bianisotropic nanoparticle 
leads to an expression known as the Sipe—Kranendonk 
conditiorF^®^. It is important to note that for a plas- 
monic bianisotropic nanoparticle, due to the intrinsic 
Ohmic losses of metals, the scattering and absorption 
cross section s will be different for froward and backward 
illuminationP^Ell, Furthermore, a planar periodic array 
of lossy bianisotropic nanoparticles, possesses interesting 
optical features such as strongly asymmetric reflectance 
and perf ect ab sorption. These effects have been stud¬ 
ied beford^^^. Although the scattering/extinction cross 
sections are similar for different illumination directions 
in the proposed nanoparticle, the contributions of the 
electric and magnetic dipole moments to the total scat¬ 
tering cross sections are not, i.e. ^ and 

[see red and blue dashed lines in Figs, ^a) and|^c)]. It 
can also be seen from Figs. ib) and[2p) that the mag¬ 
netic field distributions are significantly different when 
the nanoparticle is illuminated from froward and back¬ 
ward directions. This is obviously due to the presence 
of magneto-electric coupling, which is introduced by the 
asymmetry in the nanoparticle. 

In order to prove this, we have also calculated the 
individual polariza bility com ponents of the nanoparticle 
[see Figs. (a)-(c)The magneto-electric cou¬ 
pling aem is comparable with the electric o^ee and mag¬ 
netic amm couplings for the proposed nanoparticle [see 
Figs. 1^ (a)-(c)]. The level of this coupling can be tuned 
by changing the geometrical parameters of the nanopar¬ 
ticle, i.e. the height Hq or the diameter Dq of the par¬ 
tially drilled cylindrical air hole inside the high-index di¬ 
electric cylinder. Figure |^d) presents the magnitude of 
electric, magnetic, and magneto-electric polarizabilities 
of the nanoparticle at the maximum of magneto-electric 
coupling (|(aem|™^^) as a function of the height of cylin¬ 
drical air hole Hq. It confirms that the amplitude of 
the magneto-electric coupling can be tuned and its max- 




FIG. 3: (a)-(c) Individual polarizability components of the 
proposed nanoparticle. The geometrical parameters are: the 
height of the cylinder H = 150 nm, the height of cylindri¬ 
cal air hole Hq = 100 nm, and its diameter Dq = 150 nm. 
(d) Amplitude of the individual polarizability components at 
maximum of the magneto-electric coupling as a function of 
Hq. The maximum bianisotropy happens when the height of 
the air hole is approximately half of the height of the dielectric 
cylinder, i.e. Hq ^ — 75 nm. 


imum occurs when the height of the hole Hq is appro¬ 
priately half the height of the dielectric cylinder i7, i.e. 
Hq ^ = 15 nm. Notice, the magneto-electric coupling 

can also be tuned by changing the diameter of cylindrical 
air hole (not shown for the sake of brevity). 

As highlighted before, a lossless bianisotropic nanopar¬ 
ticle possesses identical scattering cross sections for both 
forward and backward illumination directions. Moreover, 
the forward radar cross section as shown in Fig.j^a) are 
identical for both illumination directions according to the 
definition of the normalized forward radar cross sectioiP^ 


Qrfs 


lim = 0,0^ = 0, tt) P 


r^oo A. 




47rA 

47rA 


|<^ee A em A ^me T ^ mm I 
jo^ee T ^mm | 


(6) 


However, the backward radar cross section depends on 
the illumination directions [see Fig. |^a)]. It measures 
the share of light that is directly back reflected. In or¬ 
der to show that, we start with the definition of normal¬ 
ized backward radar cross section for both illuminations 
QrbS’ backward radar cross section divided to the 
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FIG. 4: (a) Normalized backward (blue lines) and forward (red lines) radar cross sections when the nanoparticle is illuminated 
from backward (dashed lines) and forward (solid lines) directions for Hq — 75 nm. Note that the forward radar cross sections are 
identical and hence indistinguishable in the figure, (b) Radiation patterns for forward and (c) backward illumination directions 
at 445THz, respectively, (d) Normalized forward (solid lines) and backward (dashed lines) forces. — |feo|Einc|^ is the 

maximum force applied on a particle which exhibits an electric dipole moment with an incidence plane wave, (e) Schematic 
view of the different illuminations for forward and (f) backward illumination directions. 



FIG. 5: (a) A metasurface composed of an array of bianisotropic nanoparticles with Hq — 75 nm. (b) Reflection and transmission 
amplitude for two illumination directions shown in (a), (c) The same plots as in (b) but for the corresponding phases. The 
reference plane for phase is considered at z = 0 plane. 


geometrical cross sectioiP^ 

47rr^ |Esca = 0 , 6 ^ = tt, 0) 


Q 


* = lim 


RBS 


r^oo A 


|Ei„c|2 


47rA 


l^ee T ^em i ^ me 


( 7 ) 


where A = n {Dl2)‘^ is the geometrical cross section and 
D is the diameter of the nanoparticle [see Fig.[^. There¬ 
fore, for the reciprocal nanoparticle, i.e., 


the normalized backward radar cross section can 

be written aJ^ 

^RBS ~ 4:71 A ^ | • (8) 

Consequently, due to bianisotropy, the backscattering 
responses Q^gg are significantly different when the 
nanoparticle is illuminated from opposite directions [see 
Fig. B a)]. Hence, the nanoparticle possesses different 
radiation patterns for both illumination directions as de- 
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picted in Figs, [^b) and[^c). 

The dominant effect of the backscattering on optical 
scattering forces is shown for a dipolar particle without 
a magneto-electric respons^^. Since, the investigated 
bianisotropic nanoparticle shows a notable backscatterig 
difference, it is interesting to explore the optical scatter¬ 
ing forces in the realm of the bianisotropic particles. 

The optical scattering force on a dipolar particle for a 
plane wave incident field reads a^^ 



< F > 


^Re (VE*„, ■ p) + i^oRe (VHr„, ■ m) 


127reoC 


Re (p X m*) 


( 9 ) 


Plugging Eqs. § and § into Eq. § , one drives the 
optical scattering force for both illuminations: 


< {Im (age ± aem) (10) 

dzim (o^me Sz Q^mm) 

^3 ^ >1 

[(^ee Sz Q^em) (^me Sz Q^mm) ] r • 


The difference between the forces applied on the nanopar¬ 
ticle in the direction of propagation assuming reciprocal 
nanoparticle leads to 


AF = I < F+ > I - I < F- > I (11) 

= ^eo|-E’™|^Re (aemamm - "eeaL) • 

Equation 0 shows that the magneto-electric polariz¬ 
ability can lead to different optical scattering forces for 
opposite illuminations. Eigure (d) shows the normal¬ 
ized optical scattering forces. It can be seen that AF 
is pronounced when the nanoparticle exhibits a notable 
difference in backward radar cross sections (i.e. at 445 
THz). 

In summary, the following relations hold for the pro¬ 
posed lossless bianisotropic nanoparticle: 


7^ p , 7^ m , 

Csca = C'ext = C'^at = QrFS ~ ^RFS’ 

Qrbs ^ Qrbs’ I ^ > I 7 ^ I < F > |. ( 12 ) 


Now we finish the investigations of the properties of the 
individual nanoparticle and start to demonstrate its in¬ 
teresting characteristics when used in a planar periodic 
array. 


II. BIANISOTROPIC META SURFACES 

Next, we consider a periodic array composed of the 
proposed bianisotropic nanoparticle that is arranged 
along a planar surface, here the xp-plane [see Fig. da)]. 
Since the nanoparticles are assumed to be lossless, the 


FIG. 6: (a) Phase-change curve for the reflection coefficient 
of an array of nanoparticles versus the diameter of air hole at 
resonance frequency / = 405 THz for two different illumina¬ 
tion directions denoted by and r~. Ho is fixed to 75 nm. 
(b) The same plot as in (a) for variation of air hole height. 
Dq is fixed to 150 nm. 


Ohmic losses of the proposed array is also zero. Let us 
illuminate this array at normal incidence from two direc¬ 
tions, i.e. forward and backward direction [see Fig. [^a) 
for the definition of forward and backward illumination]. 
The reflection and transmission amplitudes and phases 
for these two illumination cases are plotted in Figs. |^b) 
and[^c). They were obtained from full wave simulations 
using the COMSOL MultiphysicP^. It is obvious that 
the amplitudes for the reflection and transmission for the 
different illumination directions are identical since there 
are no losses. The phase in transmission is equally iden¬ 
tical due to reciprocity. On the contrary, the phase in 
reflection is not. This is apparently due to the presence 
of bianisotropy in the proposed array. 

Now, we claim that the proposed nanoparticle is a 
proper candidate as a unit cell in a reflect arrays. In¬ 
deed, a unit cell must provide a 2^ p hase change prior 
being applicable in a reflect array^^^^^J^. This 27r phase 
change shall be obtained by a suitable geometrical tun¬ 
ing. Moreover, it has to maintain a high reflection am¬ 
plitude across the considered geometrical configuration. 
Figure shows the amplitude and phase variations ver¬ 
sus different height and diameter of the air hole inside 
the dielectric cylinder when the proposed nanoparticle is 
used in a square array with period A = 400 nm. 

As can be seen from Fig.j^a), the proposed nanoparti¬ 
cle satisfies the required condition for the phase variation 
of the reflection coefficient for both illumination direc¬ 
tions when we fix the height Hq = 75 nm and change the 
diameter = 20 — 280 nm. The amplitude of the re¬ 
flection coefficient is close to unity within an acceptable 
range while it drops down to one half at both ends of the 
diameter variations. The situation is a bit worse for the 
case when we fix the diameter Dq = 150 nm and vary the 
height Hq = 10 — 140 nm. In this case, we obtain a phase 
variation of 27r only for the forward illumination direction 
while the amplitude of the reflection coefficient, in some 
parameter regions, drops down to 20% [see Fig.j^b)]. We 
should mention that better results might be obtained if 
the geometry is carefully optimized. The presented ex¬ 
ample only serves the purpose to demonstrate the idea. 


Ho (nm) 
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Another important point is that the proposed nanopar¬ 
ticles in the array, provides asymmetric phase variations. 
That is, the variation of reflection phases are different for 
different illumination directions. It means that we may 
properly design a reflectarray with two different proper¬ 
ties when looking from different directions onto the plane. 
This is impossible using symmetrical structures which do 
not provide bianisotropic properties. 

The most important point about the proposed reflec¬ 
tarray is its transparency out of the resonance band. 
Indeed, in the investigate d reflect array, we do not rely 
on a metallic back reflectoJ ^^ l ^^ l ^^ l Instead, we have of¬ 
fered resonant nanoparticles with bianisotropic proper¬ 
ties to obtain a full reflection. This gives us the pos¬ 
sibility to preserve transparency outside the resonance 
band in the reflectarray. This transparency is very im¬ 
portant when combining multiple receiving/transmitting 
systems. Then, we need a reflectarray at a frequency 
band while we do not want to prevent signals to get trans¬ 
mitted out of that frequency band. 

III. CONCLUSIONS 

We have proposed a novel design for high-index di¬ 
electric nanoparticle which supports an omega-type bian¬ 
isotropic coupling in addition to magnetic and electric 
optically-induced dipole resonances. We have demon¬ 
strated that the magneto-electric coupling can be tuned 
by geometric parameters of the nanoparticle. In par¬ 


ticular, the nanoparticle can possess different backward 
cross sections and consequently different optical scatter¬ 
ing forces, when being illuminated by a plane wave from 
opposite directions along with the identical forward cross 
sections. 

For a metasurface created by a periodic array of bian¬ 
isotropic nanoparticles, we have observed interesting ef¬ 
fects, e.g. asymmetric reflection phases for the opposite 
illumination directions, and a possibility to achieve a 27r 
phase change together with an acceptable reflection am¬ 
plitude across the entire phase spectrum by tuning the 
geometric parameters of the nanoparticle. Finally, we 
have demonstrated that the employment of the proposed 
resonant nanoparticle together with the absence of a fully 
reflective metallic screen gives the opportunity of out-of- 
hand transparency in reflect arrays. 

All the effects described here are the direct outcome 
of the engineered omega-type bianisotropy that may 
open a new direction to design high-index nanoparti¬ 
cles and metasurfaces, including reflect arrays, transmi- 
tarrays, and Huygens metasurfaces. 
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